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Abstract 


High power rockets are essentially large unguided model rockets that can fly to altitudes 
as high as 13km and recover to earth by parachute. A rocket flight simulator can be used 
to predict the flight path of the rocket during the ascent. To do this the simulator requires 
certain input parameters. This paper summarizes some simple methods for estimating 
these parameters from very basic information such as the geometry of the rocket and the 
thrust curve. Methods for estimating the mass, centre of mass, moments of inertia, thrust 
damping, centre of pressure, drag force and normal force are given in this paper. The 
methodology described in this paper has been used in conjunction with a rocket flight 
simulator to predict the flight path of a high power rocket. Some results are given in the 


summary of the paper. 


1 Introduction 


Rocket flight simulators require certain input parameters that describe the properties of the 
rocket for the simulator. These include dynamic parameters such as the thrust, mass and 
moments of inertia of the rocket and also aerodynamic parameters such as the coefficients of 
the aerodynamic forces and the location of the centre of pressure. 

The parameters can be obtained in a number of different ways, for instance in the case 
of the aerodynamic parameters these could be obtained through wind tunnel experiments on 
the rocket or on a scale model. Alternatively they can be found by using computational fluid 
dynamics to numerically integrate the viscous and pressure forces over the surface of the rocket. 

While experimental and computational methods can provide accurate estimates of the pa- 
rameters it is also useful to be able to estimate the parameters using simple equations derived 
from theoretical relationships using some simplifying assumptions. This approach provides a 
quick method for estimating parameters in applications where the most accurate estimates are 
not required. 

This document collates and summarizes relevant equations for estimating the dynamic and 
aerodynamic input parameters from rocketry and engineering literature. In some cases equa- 
tions are semi-empirical, that is they use relationships that have been fitted to data from 


experiments using rockets of a standard shape. In general, these equations are only applicable 


to rockets with a conventional shape. That is an axisymmetric rocket with a single cylindrical 
body, a conical, parabolic or ogive shaped nose cone and fin sets consisting of three or four 
trapezoidal fins. Changes is body diameter via conical transitions are allowed. 

There are also a number of simplifying assumptions that are used in the estimation of 
the parameters and in the rocket simulators that use these parameters. These include the 
assumption that the rocket is an axisymmetric rigid body and that the angle of attack during 
the flight is always low. These and other assumptions will be discussed in more detail throughout 
the paper. 

The methods presented here have been used to provide estimates of input parameters for a 
rocket flight simulator described in Box et al [2009], which also contains some validation data 


from a test flight. 


2 The Dynamic Parameters 


The dynamic parameters of the rocket that need to be estimated for flight simulations are the 
Thrust, mass, centre of mass, moments of inertia about the three principal axis and a thrust 
damping coefficient. Because the rocket is expelling mass from the motor nozzle all of these 
parameters will be time-varying until the motor burns out. 

The thrust curve of the rocket motor can be measured on a test rig but in most cases 
commercially sold motors will be supplied with thrust-time data from the manufacturer. If 
the mass of the fuel is known then it is possible to estimate the mass curve for the motor by 


assuming that the mass burned at any point in time is proportional to the impulse of the motor 


m; f T dt 
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where My is the total mass of fuel, T is the thrust and ôM; is the change in the mass of fuel 


up to that point (1). 


ôM; = — (1) 


on board from ignition to time t = i. This calculation can be done numerically. 

The following sections show how to estimate the centre of mass, moments of inertia and 
thrust damping coefficient for a given rocket mass M, which constitutes a single point on the 
mass curve. So to apply these calculations they must be repeated over the range of masses to 


build up time varying curves for each of these parameters also. 


2.1 Centre of Mass and Moments of Inertia 


For the purposes of these calculations it is assumed that the rocket is an axisymmetric body 
and the centre of mass of the rocket and it’s component parts lies on the axis of symmetry, 
which is the roll axis of the rocket. We define the position of the centre of mass as a distance 
Xem along the roll axis measured from the tip of the nose cone. The position of the centre of 
mass is defined using (2). 

dm (2) 


where r is the distance of an infinitesimal “slice” of rocket with mass dm along the roll axis 
and M is the total mass of the rocket. 


Similarly the moment of inertia of the rocket about an axis of rotation is defined as 


T= [Pam (3) 


where r is the distance of an infinitesimal volume with mass dm from the axis of rotation. 

A simple approach to solving (2) and (3) is to divide the rocket up into component parts. 
The parts are then assumed to have simple shapes such as tubes, cones, blocks, cylinders etc. 
Small parts can be modelled as a point mass. Parts with complex shapes such as the packed 
parachute can approximated by a simpler shape, in this case a cylinder with uniform density. 

The following subsections describe how to calculate the centre of mass and moments of 
inertia about the principle axes for some simple shapes and how to combine them together to 


estimate the total centre of mass and moments of inertia for the rocket. 


2.1.1 Cone 


The centre of mass of a solid cone is located on the axis of symmetry two thirds of the distance 
between the tip and the base. 
2l 
Xeom= a (4) 
where Xem is distance from the cone’s point to it’s centre of mass and / is the height of the 
cone. The moments or inertia about the principal axes of a cone where the z axis is aligned 


with the cones axis of symmetry are given by 
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where M is the mass of the cone and r is the radius of the base of the cone. 


2.1.2 Tube 


The centre of mass of a tube is located half way along the axis of symmetry. 


l 
Kem = 
l (7) 


where Xem is the distance of the centre of mass from one end of the tube and 1 is the length of 
the tube. The moments of inertia about the tube’s principal axes, where the z axis is the axis 


of symmetry are 


Iss = Lyy = aM [3 (ri +73) +27] (8) 
1 
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where M is the mass of the tube and rı and ro are the internal and external radii of the tube 


respectively 


Block 


The centre of mass of a block with sides of length X, Y and Z will be at the centre of the block 
(X/2, Y/2, 2/2). 


The moments of inertia about a block’s principal axis are given by 


1 

Ing = Mv? + 2°) (10) 
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where M is the mass of the block. 


Point Mass 


Centre of mass of a point mass is it’s location. A point mass has zero moment of inertia about 


it’s own centre of mass. 


Total Centre of Mass and Moments of Inertia 


The location of the overall centre of mass of the rocket with reference to the nose cone tip is 


calculated using the following equation. 


5 (Xp zF Xem(P))Mp 
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Xem(R) = (13) 


Mr 


where the subscript P denotes a part and R denotes the whole rocket. Xp is the distance from 
the nose tip to the most forward point of part P. 


The overall moments of inertia are calculated using the parallel axis theorem. 


Ip = Temcp) + Mpd? (14) 
In = Ip (15) 
PER 


where Ir is the moment of inertia of the rocket about a given axis through the rocket’s centre 
of mass, Ip is the component of the rockets moment of inertia about the selected axis that is 
due to a given part P. Iem(p) is the moment of inertia about a parallel axis through the centre 
of mass of the part P, d is the distance between the two axes and Mp is the mass of the part. 

Note that although a point mass has no intrinsic moment of inertia about it’s own centre 
of mass it may still contribute to the moment of inertia of the rocket through the parallel axis 


theorem. 


2.2 Thrust Damping 


The expulsion of mass from the rocket nozzle has a further effect on the rocket’s dynamics and 


this is to damp angular rotations of the rocket by generating a torque in the opposite direction 


to the rotation. This is caused by the lateral acceleration of the hot rocket gas as it travels 
down the motor tube [Duncan & Ensey, 1964]. 
The torque due to thrust damping is given by 


Taa = —M (l, — 2w (16) 


where M is the mass expulsion rate of the fuel, len is the distance between the centre of mass of 
the rocket and the nozzle exit, lec is the distance between the centre of mass of the rocket and 


the centre of mass of the fuel and w is the angular velocity of the rocket. The mass expulsion 


dM 
dt |t= 


The term M(I2,, — 12) is a dynamic parameter that must be passed to the simulator. We 


rate at time t = x is , and can be estimated numerically from the mass time curve. 


can define this term as the trust damping coefficient (Caa). 


Caa = M(Èn — lao) (17) 
Tda = —Caqw (18) 


However it should be noted that while it is conventional to call this a coefficient it is not 


2 
dimensionless, it has the SI units kau 


3 The Aerodynamic Parameters 


3.1 Force Components 


Figure 1: Showing the vectors of drag force Fp and lift force Fz relative to a rocket’s apparent 
velocity vector V for flight with angle of attack a. 


The aerodynamic force vector on a rocket in flight is typically considered as two orthogonal 
force vector components, however there are two conventions for how to orientate these com- 
ponents. Figure 1 shows the first convention. Here the components of aerodynamic force are 
the drag force Fp and the lift force Fg. These aerodynamic forces can be assumed to act on a 


single point in the rocket known as the centre of pressure. The vectors of the component forces 


Figure 2: Showing the vectors of axial force F4 and normal force Fy relative to a rocket’s roll 


axis. 


are defined relative to the rocket’s apparent velocity vector V. The apparent velocity is defined 
as velocity of the centre of pressure of the rocket relative to the atmosphere. The direction of 
the the drag force vector is exactly opposite to the direction of V. 

The second convention is to orient the two orthogonal components relative to the rocket’s 
roll axis as shown in Figure 2. Here the forces are referred to as the axial force F4 and the 
normal force Fy. As with drag and lift these forces act through the centre of pressure. 

The angle between the vector defined by the rocket’s roll axis and the apparent velocity 
vector V is known as the angle of attack a. The direction in which lift force or the normal force 
acts is orthogonal to the drag force or the axial force and in the plane defined by the roll axis 
vector and the apparent velocity vector. Also in both conventions the normal force, or the lift 
force goes to zero as the angle of attack a goes to zero. 


The magnitude of any of these forces can be given by the following equation 


1 


where the subscript i is either D,L,A,N for drag, lift, axial or normal respectively. p is the 
atmospheric density V is the rocket’s apparent velocity A, is the reference area, which is 
the cross-section area at the base of the nose cone and C; is the relevant aerodynamic force 
coefficient. The challenge is to estimate these aerodynamic coefficients. 

The aerodynamic forces are complex phenomena and there are no simple analytical solutions 
for the coefficients. In the incompressible flow regime the forces can be divided into pressure 
force and viscous force. Pressure force arises through the stagnation of fluid on the rocket 
forebody, fins, and any other protrusions, and also through a suction force created by a low 
pressure region at the base of the rocket where boundary layer separation occurs. Viscous force 
is due to skin friction between the rocket and the air. If the boundary layer at the rocket’s 
surface is turbulent then the viscous forces will be significantly less that if it is laminar, hence 


the viscous force is highly dependant on Reynolds number. 


As the rocket approaches the speed of sound a further force becomes relevant. At trans- 
sonic speeds shock waves form at the nose tip and at the leading edge of the fins. At sonic 
velocity these shocks are normal to the direction of motion of the rocket, and they become 
oblique (decreasing in half angle) as the rocket proceeds through super-sonic speeds. Momentum 
is transferred from the rocket to the surrounding air via these shock waves. This increased 
component of aerodynamic force is sometimes called wave drag. 

At low speed (incompressible flow) the aerodynamic coefficients are functions of the angle of 
attack (a) and Reynolds number(Re). At higher speeds where the flow becomes compressible 
(Ma > 0.4) they are also functions of Mach number (Ma). 


Ci = f(a, Re, Ma) (20) 


3.2 Equations for estimating the Aerodynamic Coefficients 


This section presents two separate methods from rocketry literature for estimating an aero- 
dynamic coefficient using the rocket’s geometry. The first is for estimating the normal force 
coefficient Cy and the second is for estimating the coefficient of drag Cp. 

It might be considered more helpful to present two methods based on the same convention 
however this is not available in the literature. It is of course possible to combine these two 
methods to make estimates of either Cy and C4 or Cp and Cz. This however leaves the 
problem that two separate methods published by different authors are being used to estimate 
two components of what is essentially a single force, the aerodynamic force on the rocket. 

Despite these misgivings, test data from real rocket flights has been shown to compare well 
with flight simulation data based on this method (see Section 4). 

Figure 3 shows a schematic of an example rocket shape. On this figure key dimensions that 
are used to define the rocket’s geometry are labelled. This figure will be referred to throughout 
the following sections. 


3.2.1 Normal Force 


Equations for estimating the normal force, and the location of the centre of pressure on a small 
rocket have been proposed by Barrowman [1998]. The Barrowman equations are derived using 


the following assumptions. 
e The angle of attack of the rocket is low (a < 10°) 
e The effects of compressibility can be neglected i.e. Ma < 0.4 
e Viscous forces are negligible 
e Lift forces on the rocket body tube can be neglected 
e The air flow over the rocket is smooth and does not change rapidly 
e The rocket is thin compared to it’s length 


e The nose of the rocket comes smoothly to a point 
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Figure 3: Schematic of a simple rocket with showing the principle dimensions 


e The rocket is an axisymmetric rigid body 
e The fins are thin flat plates 


Here we will summarize the Barrowman method, a full derivation is given in Barrowman & 


Barrowman [1966]. From (19) the normal force on the rocket is given by (21) 


1 
Fy = eV ArCn (21) 


where A,. is the cross sectional area of the rocket at the base of the nosecone. The Barrowman 
method assumes incompressible flow and neglects viscous forces, therefore Cy is assumed to 
be a function of a only. Furthermore because a small œ is assumed Cy can be expressed as a 


linear function of a (22). 


Cn = Cnoa (22) 


where C'yq is the slope of the normal force coefficient, sometimes called the stability derivative. 

CNa is calculated as the sum of the individual stability derivatives for each of the components 
that makes up the rocket’s shape. These components can be either a nosecone, a finset or a 
conical change in the diameter of the rocket’s body tube. The body tube itself is neglected 
as a component because one of the key assumptions of the Barrowman method is that lift on 
the body tube is negligible. A rocket can consist of a nosecone and any number of finsets and 
conical changes in body diameter. 

The total stability derivative is given by the sum of the individual stability derivatives for 
each of the parts of the rocket. 


CNa(R) = 5 CNa(P) (23) 
PER 


where P designates a rocket part and R designates the whole rocket. 

The location of the rocket’s centre of pressure is also calculated is terms of the centre of 
pressure of each of the individual rocket components. As with the centre of mass it is assumed 
that the centre of pressure lies on the roll axis of the rocket, the position of the centre of 
pressure is defined as the distance (Xep) of the centre of pressure from the nose cone tip. Xep 
is calculated using equation (24) 


XO CnaP) Xepe) 
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where Xep(p) is the distance of the centre of pressure of component P from the tip of the 
rocket’s nosecone. 
The equations for calculating Cyna and Xep for nosecones, finsets and conical changes in 


body diameter are given below. 


Nose cone The stability derivative for a rocket’s nosecone is 2 as long as the shape of the 


nosecone is either conical, ogive or parabolic. 


C'Na(n) =2 (25) 


The location of the centre of pressure of the nose cone is dependent on the shape 


2 
Conical: Xep(n) = gin (26) 
Ogive: Xep(n) = 0.4661, (27) 
1 
Parabolic: Xep(n) = gin (28) 


where Xep(n) is the distance of the centre of pressure from the tip of the nose cone and I, is 


the length of the nose cone. 


Conical change in body diameter The stability derivative Cya and the location of the 


centre of pressure Xep for a conical change in body diameter are given by the following equations. 
Ga (daN? 
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Xoe) = Xet $ |1+ a (30) 
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where dn is the diameter of the base of the nose cone (reference diameter), du is the upstream 
diameter of the conical change, dg is the downstream diameter of the conical change, le is the 
length of the conical change and Xe is the distance between the tip of the nose cone and the 

most upstream point on the conical change in body diameter (figure 3). 
Equation (29) will produce a negative result if the upstream diameter is larger than the 
downstream diameter. This is correct as the pressure force in this case is a suction force acting 


towards the direction of flow. 


Fins The following equations are for the entire fin set, and are valid for configurations with 
3 or 4 trapezoidal fins. 
The stability derivative for a trapezoidal finset is given by the following equation 


2 
4n (4) 
ae ee: (31) 
1+ 1+ (24) 


where n is the number of fins, ls is the fin span, lm is the fin mid-chord, l, is the fin root-chord 


Cna(f) = K fo 


and l; is the fin tip-chord (figure 3). Arp is a coefficient that takes into account an increase in 
the normal force due to interference effects between the fin and the body. 


d 
Kp = 1+ — (32) 
(ls + £) 


where dy is the diameter of the body tube at the fin’s location. The location of the centre of 


pressure for a trapezoidal finset is given by 


1 Lol 
Xop = Xpt BEET) Efren- a (33) 


where Xy is the length between the nose cone tip and the point where the fin leading edge 
meets the body tube. 


Expanding the Valid Range of the Barrowman Equations As mentioned in the above 
section the Barrowman equations are subject to some quite restrictive assumptions. At least 
two of these assumptions may be broken in high powered rocket flight. The first of these is 
the assumption of incompressibility as high power rockets are capable of exceeding the speed 
of sound. 

The second assumption is that the lift forces on the rocket body are negligible because 
the angle of attack a is below 10°. This assumption may be reasonable for most rockets but 
experiments have suggested that for rockets with particularly long slender bodies, body lift may 
not be negligible [Dahlquist, 1998]. 

Correcting for the effects of compressible flow is dealt with in section 3.2.3. Below we present 


a method for estimating the lift on the rocket body. 


Rocket body lift correction Rocket body lift and it’s associated effects on the position of 
a rockets’ centre of pressure are neglected in the Barrowman equations. Galejs [1999] proposed 
the following extension to the Barrowman equations to take body lift into account. 

The coefficient of normal force due to body lift is given by equation (34). 

Cnu) = KPa? (34) 
r 

where K is a constant between 1.0 and 1.5, A, is the planform area of the rocket (excluding 
the fins) and A, is the reference area for the rocket (cross-sectional area at the base of the nose 
cone). 

C'y(z) is not a linear function of a as with the assumption of the Barrowman method. Galejs 
defines 


A 
Cno2 = Ke (35) 


CN(L) 
Q 


which is not the slope of Cy) but simply . CNa2 can be added to Cya in equation (23). 
But this makes Cya a function of a and hence the whole rocket Cy is non-linear. 

The centre of pressure for the body lift force is taken to be the centre of the planform area, 
The equations used to determine the centre of the planform area for the various rocket body 


components are given below. 


Conical Nose: Xep(n) = In (36) 
Ogive Nose: Xep(n) = ty (37) 
Parabolic Nose: Xep(n) = ly (38) 
Body section: Xep(n) = Xe + sl (39) 


Conical transition: X, 


aches z (40) 


= Xe — | ———_ 
cp(n) + 3 | (da + da) 
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For each component Cy,2 and Xep are included in the Barrowman equations (23) and (24). 
In the Galejs [1999] article results of this modelling were compared with experimental data 
for variation in CP with angle of attack from Dahlquist [1998]. The comparison shows a good 


agreement between experimental and modelled data for values 0° < a < 15°, using K = 1. 


3.2.2 Drag Force 


Equations for estimating the coefficient of drag force on a rocket are provided by Mandell et al 
[1973]. First we present how to calculate the drag force at zero angle of attack, then we show 


how so extend this to small angles of attack using semi-empirical equations. 


Zero angle of attack drag force The coefficient of zero angle of attack drag force (C'p(o)) 
on a rocket can be estimated using the United States Air force Stability and Control Datcom 
Method [Mandell et al, 1973]. This method, like the Barrowman equations, uses basic informa- 
tion about the rocket geometry and assumes incompressible flow. to estimate the drag forces. 
The main equations of the datcom method are summarised below; derivations of the equations 
are given in Mandell et al [1973]. 


Body Drag The drag on the rocket forebody is estimated using equation (41) 


60 ly ln ly da le 
= |1 + ——_ .0025—| |2.7— + 4— +2 [{ 1- — |] =I C 41 
CDi) T (IrrR/d)’ ii =] | dy = dy : ( z) a PLO) oe 


where [7p is the total length of the rocket body, le is the length of a boat tail (if present) - a 
boat tail is a conical reduction in the body diameter at the base of the rocket (figure 3). dẹ is 
the maximum body diameter and d4 is the diameter of the rocket base. Cf sy) is the coefficient 
of viscous friction on the rocket forebody (defined later in (45)) 


Base Drag The base drag on the rocket is the drag due to the low pressure region at the 
base of the rocket that is caused by boundary layer separation. This drag is estimated using 
equation (42) 


(2) 


Cp) = 0.029— (42) 
VCD(fo) 
Fin drag The fin drag on the rocket at zero angle of attack is given by equation (43). 
Ts \ 4nA 
Cory = 2C pp (14274 ) SP (43) 
in td F 


where Cpp) is the coefficient of viscous friction on the fins (defined later in (45)), Ty is the 
fin thickness, n is the number of fins and dy is the diameter of the body tube at the fin root. 

Afp is defined as the fin planform area. The planform area of the exposed part of a trape- 
zoidal fin is given by Afe = 4 (Ip +1,)l,. This is known as the exposed area. For the full 
planform area it is usual to assume that the fin extends to the centre line of the rocket body 
ie. Afp = Afe + $dflr. 
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Interference Drag The drag due to interference effects between the fins and the body is 
given by equation (44) 


Tr\ 4n(Ar,y — Afe 
Cpi) = 2C Ff) (1 + 272) an (44) 


Viscous Friction As discussed in section 3.1 the viscous forces on the rocket are dependent 
on Reynolds number. The friction force coefficient is given by equation (45). 


= when Re < Re, 
Cy = (45) 
rare — = when Re > Re, 
where B is given by equation (46). 
B = Reg (San = =| (46) 
Reynolds number is given by 
Re= ane (47) 


where p is the atmospheric density, js is the kinematic viscosity of air, V is the apparent 
velocity vector and L is the characteristic dimension. To calculate the viscous friction coefficient 
of the rocket body C(») the characteristic dimension is the total body length L = Irr. To 
calculate the viscous friction coefficient of the rocket fins C's) the characteristic dimension is 
the fin mid chord L = lm. Re, is the critical Reynold number, which Mandell et al [1973] give 
as 5 x 105. 


Total zero angle of attack drag The total zero angle of attack drag coefficient Cp(o) can 
be modelled as sum of the individual drag coefficients given in equations (41) to (44) 


Cpo) = Cafe) + Cpo) + Cpi) + Cac) (48) 


Additional drag at angle of attack Mandell et al [1973] presents equations for estimating 
an additional component of drag, which can be added to the zero angle of attack drag (C'p(o)) 
in order to model the drag force at small angles of attack a. This method incorporates some 
coefficients that are derived from data which comes from wind tunnel experiments on rocket 
models. 

As with the zero angle of attack drag force the component of alpha drag is subdivided into 
drag on the rocket body Cppq) and drag on the finset Cp f(a). 

The coefficient of alpha drag on the rocket body is calculated using equation (49). 


3.67(1.36lrR — 0.551, 
Cpo(a) = 25a? + Seen USBI) (49) 
Tdp 
where a is the angle of attack of the rocket and 6 and ņ are experimentally derived coefficients. 
These are obtained using experimental data from wind tunnel measurements. Figure 4 from 


Mandell et al [1973] shows ô and n 
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Figure 4: Graphs showing the relationship between a and the coefficients 6 and 7 [Mandell et 
al, 1973]. 


The coefficient of alpha drag on the rocket’s fins is calculated using 


ded Ajed 
: — 1 
Td? +3.12(k fo + kof (2 


CDfla = 0° |1.2 (50) 


where Afp is the total fin planform area and Afe is the total fin exposed area. krp and kps are 
the fin-body interference coefficient and the body-fin interference coefficient given by equations 
(51) and (52). 


II 


kyo 0.8065R? + 1.1553R, (51) 
kop = 0.1935R2 +0.8174R, +1 (52) 


where Rs is the fin section ratio, which is the ratio between the total span on the fins lrs (figure 
3) and the diameter of the body tube where the fins are mounted dy, (Rs = FS). 
The complete coefficient of drag on the rocket is obtained by adding the alpha drag coeffi- 


cients to the zero angle of attack drag coefficient, hence 
Cp = Cpo) + Cdbva + CDfa (53) 


Knowing estimates for the drag force coefficient Cp and the normal force coefficient Cy an 
estimate for the axial force coefficient C'A can be calculated using 
Cp cosa — £Cn sin (2a 
ae (54) 
1—sin* a 


This gives us two orthogonal force coefficients using the same convention (Cy and C4). 


3.2.3 Compressible Flow Correction 


The methods described above for calculating Cy and C'4 are only valid for incompressible flow. 


In this section we show how to extend the range of these aerodynamic coefficients into the 
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compressible flow regime. The Prandtl-Glauert compressibility correction can be applied to 
estimate the variation in the aerodynamic force coefficient at trans-sonic speeds. 
At subsonic speeds (Ma < 1) the corrected aerodynamic coefficient is given by equation 
(55) (Cramer, 2002] 
gds- (55) 
V1 — Ma? 
where C; is the aerodynamic coefficient in the incompressible regime and Ma is the free stream 
Mach number. 
At supersonic speeds (Ma > 1) the corrected aerodynamic coefficient is given by equation 
(56). 
deN (56) 
VMa?—1 
The problem with this estimation method is that C? goes to infinity as Ma goes to 1, which 
no longer fits experimental data. It is suggested in Ketchledge [1993] that in order to avoid this 
situation, the following equation is used for 0.8 < Ma < 1.1. 


Io Ci 
= i= (0.8)2 a 


4 Summary 


In this document it has been shown how to calculate a number of dynamic and aerodynamic 
parameters for rocket flight simulators. The parameters which have been dealt with are sum- 
marized in Table 1. As the table shows, many of the parameters can be functions of a variable 
in the rocket simulation such as time t, Reynolds number Re or angle of attack a. There are 
different ways that these dependencies can be dealt with, the parameters can be calculated on 
the fly for the simulator, or alternatively pre-calculated at discrete intervals and stored in a 


database to be accessed by the simulator and interpolated as required. 


Thrust 
Mass 


Distance of the centre of mass from the nose tip 


Moments of inertia about the rocket’s yaw axis 


Moments of inertia about the rocket’s pitch axis 
Moments of inertia about the rocket’s roll axis 
Thrust damping coefficient 


Coefficient of normal force 


Re, a, Ma Coefficient of aerodynamic roll torque 


Distance of the centre of pressure from the nose tip 


Table 1: Parameters that are stored in the database for the rocket model. 
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The list of parameters that has been dealt with is not exhaustive. These parameters provide 
enough information for a simulator which assumes that the rocket is an axisymmetric rigid body. 
However for more complex simulations more parameters may be required, for example when 
simulating a non-axisymmetric rocket or modelling other phenomena that are neglected here 
such as thrust misalignment. 

Nevertheless simulations using the axisymmetric assumption and the parameters addressed 
in this paper have been shown to be useful. Box et al [2009] used the methods laid out in this 
paper to provide inputs to a rocket flight simulator which was used to simulate the flight of a 
high powered rocket. The rocket was flown to an altitude of around 3.5 km before descending 
under parachute. Instrumentation on board the rocket recorded some data on the flight and 
the landing position was logged with a GPS receiver. Table 2 and Figure 5 from Box et al 
[2009] show some comparison between simulated and measured data. Various flight statistics 
are given in Table 2 and Figure 5 shows the flight path of the rocket ascent and parachute 
descent as predicted by the simulator. The ellipses in the figure mark the areas of lo and 2c 
confidence in landing position. The diamond marker shows the recorded landing position of the 
actual rocket. 

This is a single result and not a comprehensive test of the simulator or the methodology 
used to estimate the input parameters, but it is an encouraging result that suggests that future 


investigation into the validity of this methodology may be worthwhile. 


COo o i o a 


Landing position |E, N] [-135m, 936m] | [—71m, 1042m] 


Difference in landing positions 


Table 2: Comparison between simulated and measured flight statistics. 
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Figure 5: Mean simulated flight path with lo and 20 landing probabilities. The measured 


landing position is marked by a diamond. 
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